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Abstract 

Zernike polynomials are commonly used to represent the wavefront phase on circular 
optical apertures, since they form a complete and orthonormal basis on the unit disk. In 
[Diaz et all, 2014] we introduced a new Zernike basis for elliptic and annular optical aper¬ 
tures based on an appropriate diffeomorphism between the unit disk and the ellipse and the 
annulus. Here, we present a generalization of this Zernike basis for a variety of important 
optical apertures, paying special attention to polygons and the polygonal facets present in 
segmented mirror telescopes. On the contrary to ad hoc solutions, most of them based on 
the Gram-Smith orthonormalization method, here we consider a piece-wise diffeomorphism 
that transforms the unit disk into the polygon under consideration. We use this mapping to 
define a Zernike-like orthonormal system over the polygon. We also consider ensembles of 
polygonal facets that are essential in the design of segmented mirror telescopes. This gener¬ 
alization, based on in-plane warping of the basis functions, provides a unique solution, and 
what is more important, it guarantees a reasonable level of invariance of the mathematical 
properties and the physical meaning of the initial basis functions. Both, the general form 
and the explicit expressions for a typical example of telescope optical aperture are provided. 

2010 AMS Mathematics Subject Classification: 41A10 (Approximation by polynomials); 41A63 
(Multidimensional problems); 78M99 (None of the above, but in this section: Optics, electro¬ 
magnetic theory); 85-08 (Computational methods: Astronomy and astrophysics). 

Keywords & Phrases: Zernike polynomials; orthonormal systems; polygonal facets; segmented 
mirror telescopes. 


1 Introduction 

The problem of finding complete orthononnal systems to represent functions defined on finite 
supports with a given geometry appears in many areas of Physics and Engineering. In particular, 
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Zernike circle polynomials [31] are widely used to represent optical path differences (phase 
differences or wave aberrations) in wavefronts, or even the sag of optical surfaces (such as the 
human cornea [28], [9]) as they are well adapted to the circular shape (disks or slightly deformed 
disks) of a majority of conventional optical systems. There is an infinite number of possible 
systems, but Zernike polynomials (ZPs) (or lineal combinations of them [11]) show important 
advantages and interesting properties. Among these properties, ZPs permit to establish a link 
with the traditional Seidel theory of aberrations [1], which is based on a third order Taylor series 
expansion, and with further extensions of the Seidel theory to 5th order, etc. The theoretical 
and practical advantages of orthonormal polynomials, make that ZPs became the standard way 
to describe the phase of wavefronts [22] (or the wave aberration or optical path differences) in 
many fields ranging from atmospheric optics [23], optical design and testing [17] or visual optics 
(the ANSI Z80.28 standard for reporting aberrations in the human eye is based on ZPs). 

The circle is the most common optical aperture, but other more complicated geometries 
have gained interest in recent years after the development of large telescopes that require more 
sophisticated optical designs such as the annular pupils in large telescopes [26]. In particular, 
segmented mirror telescopes (SMTs) are commonly used nowadays by NASA, ESA and other 
astronomical and astrophysical organizations all over the world to collect information from the 
outer space in the form of electromagnetic radiation. An SMT is an array of smaller mirrors 
designed to act as segments (or facets) of a single large mirror. The facets can have diverse 
shapes, ranging from planar squares to curved asymmetric polygons [24]. They are used as 
objectives for large reflecting telescopes. To function, all the mirror segments have to be polished 
to a precise shape and actively aligned by a computer controlled active optics system using 
actuators built into the mirror support cell. 

An important example of SMT is The European Extremely Large Telescope [5]. Because 
current monolithic mirrors cannot be constructed larger than about eight meters in diameter, 
the use of segmented mirrors is a key component of current large-aperture telescopes. This is 
because of the technological limit of a primary mirror made of a single rigid piece of glass [21]. 
Using a monolithic mirror much larger than 5 meters is prohibitively expensive due to the cost 
of both the mirror, and the massive structure needed to support it. A mirror beyond that size 
would also sag slightly under its own weight as the telescope was rotated to different positions 
changing the precision shape of the surface [7]. Facets are also easier to fabricate, transport, 
install, and maintain over very large monolithic mirrors. 

Other examples of SMTs are: (i) The Keck Telescope in Hawaii, the largest optical and 
infrared telescope in the world, that consists of 36 hexagonal segments, each 0.9 meters on 
aside (http://www.ps.uci.edu/physics/news/chanan.html). (ii) The Gran Telescopio Canarias 
(GTC), whose useful wavelength range extends from 365nm to 25 /jm, its primary mirror consists 
of 36 hexagonal facets of Zerodur coated with aluminium; each segment measures 1.9m from 
vertex to vertex and has a side length of 0.9m and a weight of 470 kg. They cover a total 
surface of 75.7 square metres and have gaps of 3mm between them (http://www.gtc.iac.es/ 
observing/GTCoptics,php$\#$Ml). 

The importance of these non-circular geometries motivated the development of a series of 
ad hoc solutions, most of them based on the Gram-Smith (G-S) method to obtain orthonormal 
basis on different types of apertures [29], [30] such as ellipses [4], rectangles [15], annuli [12], 
circular sectors [3], etc. The main drawback of the G-S method is that it does not provide of a 
unique solution that hinder the physical interpretation of the associated expansion coefficients 
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Figure 1: Example image of a wavefront test (left) and example of wave- 

front solution (right). Pictures taken from the GTC theoretical base data: 
http: / / www.gtc.iac.es / observing/GTCoptics.php#Ml. 


(especially for the higher orders due to a cumulative effect associated to the G-S method). 

In our previous paper [20], we proposed a mapping that transforms the unit circle into the 
desired geometry, paying special attention to elliptic and annular aperture geometries. Using 
this mapping, we obtained a new Zernike basis for elliptic and annular optical apertures with 
important mathematical and physical properties. In this paper we go a step further giving a 
rigorous formulation under a unique criterion and providing a unique general solution for most 
of the usual optical apertures. As it is explained in [20], our approach is based in finding an 
appropriate mapping that transforms the unit circle into the desired aperture geometry. That 
is finding piece-wise diffeomorphism that transforms the unit circle into the set within the plane 
that represents the optical aperture. This mapping means warping the input basis functions so 
that they fit into the new aperture geometry. On the contrary that ad hoc solutions, that warping 
permits not only unicity, but also a high level of invariance of the mathematical properties and 
physical meaning of the basis functions (tilt, defocus, astigmatism, coma, etc.) and hence a 
natural generalization of the aberration theory. 

Although this technique may be applied to different optical systems defined by different 
geometries (ellipses, circular sectors, annuli,...), in this paper we focus our attention in the 
geometries required for the SMTs: polygons and ensembles of polygonal facets. Apart from 
SMTs, a different optical context where this kind of geometry is of interest can be found in the 
high resolution compound superposition eyes of some diurnal insects that suggest high resolution 
optics [19]. 

In the following section we report the definition of the Zernike orthonormal system and 
revise, within a rigorous formulation, the general method introduced in [20], that is valid for 
any domain, piece-wise diffeomorphic to the unit disk. In Section 3 we focus our attention on 
polygonal domains and give a detailed description of the corresponding Zernike-like system. In 
Section 4 we consider an ensemble of hexagons typically used in the primary phase of SMTs and 
give an orthonormal system composed by a combination of the Zernike-like systems of the facets. 
In Section 5 we give an example of wave front description in a polygon and in an ensemble of 
polygonal facets using the Zernike-like basis introduced in Sections 3 and 4. Some final remarks 
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are given in Section 6. 


2 Zernike-like systems in compact sets of the plane 


2.1 Zernike polynomials in the unit disk 


Through out this paper, u := (u,v) represent cartesian coordinates in the unit Disk D and p, 4> 
the corresponding polar coordinates: (u,v) = a(p,(j>) := (pcoscj), psincj)), 0 < p 2 = u 2 + v 2 < 1 
and 0 < 4> = arctan (y/x) < 2n. Also, we will use the notation Z™(u,v ) for the Zernike 
polynomials in cartesian coordinates and the notation Z™(p, (p) for the Zernike polynomials in 
polar coordinates defined in the unit disk D: Z™ := Z™ o a. The standard definition of the 
Zernike polynomials in polar coordinates is the following [22]: 


z™(pA) 
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— N™R\™\p) sin(m0), m < 0. 
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where 5 m p is the Kronecker delta function: 5 q,o = 1 and o = 0 for m / 0. The index n 
indicates the degree of the radial polynomial R^(p) and the index m the azimutal frequency of 
the azimutal component. The index n is any positive integer or zero: n = 0,1, 2,.... For a given 
n, the index m takes the values m = — n, — n + 2, — n + 4, ..., n. More details, as the expression 
of the Zernike polynomials in cartesian coordinates may be found in [20] for example. For later 
convenience, we summarize here the orthogonality and completeness properties of the Zernike 
polynomials: 

The family { Z ™} n + 2 ’—’" j s an orthonormal system in L 2 (D). Here, and in the re¬ 
maining of the paper, L 2 (D) represents the set of square integrable functions defined on the 
unit disk D with respect to the normalized Lebesgue measure dp := dudv/ir. Equivalently, 
the family n+2 ’"'’ Tl is an orthonormal system in L 2 (D ) with respect to the measure 

dp := pdpdcf)/TT [ [22], chap. 2]: 


/ rym ryr 


D 


Z™.Z%'dp = (Z™,Z%') 




z% dp = 5, 


n,n' 
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Moreover, the Zernike set is a complete orthomormal system of L 2 (D). That is, any square 
integrable function / : D —>• C can be written in the form [ [22], chap. 2]: 

oo n oo n 

/ = £ E c K, f = T,T, (3) 

n=0m=—n n=0 m=—n 

where / := / o a is the function / in polar coordinates. In the above equalities: 


C:= (Z^J)^ 


Z™-fdp, 
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The equality (3) must be understood in the L 2 sense, that is, the convergence of the series in 
the right hand side is understood in the L 2 (D) norm ||/|| 2 := (/, /) /i = (/, f)p,: 


lim 

TV —yoo 


N n 

/-EE « 

n= 0 m=—n 


lim 

N —>oo 


N n 

n=0 m=—n 


= 0 . 


2.2 Zernike-like systems in other sets of the plane 

Any planar optical aperture may be mathematically described as a set M in the plane piece-wise 
diffeomorphic to the unit disk D. Then, we consider all the subsets M C M 2 obtained from D 
by means of a mapping p : D M that is a diffeomorphism a.e. in D (see Figure 3): 

ip : D —> M CM 2 , 

(u,v) ->• (x,y) = p(u,v). 



Figure 2: The unit disk D and the set M are deformed each in other by the respective mappings 
p and p~ x . 


In the remaining of the paper we use the notation x := (x,y) for the cartesian coordinates 
in M and r 2 = x 2 + y 2 and 6 = arctan(y/x) for the corresponding polar coordinates in M, with 
0 < 8 < 2-k. Then, for (x,y) € M we write x := (x,y) := p(u,v) = (x(u,v),y(u,v)). And 
reciprocally, u := (u,v) = p~ 1 (x,y) = (u(x,y),v(x,y)). We denote by J(x,y) the jacobian of 
</? -1 that is defined a.e in M: 


J(x,y) 



and dpj := \J(x,y)\dxdy/'K the measure induced by p in M. In the remaining of the paper we 
assume that J(x,y) is continuous a. e. in M. Then, it is clear that: 




dp = 


(z; 


m op- 1 ) • (Z™ op- 1 ) dpj. 


J JD J JM 

In other words, the set of functions 

AT™ := Z™ o p- 1 ; n = 0,1,2,...; m = — n, — n + 2, — n + 4,... , n, 

defined in M, is an orthonormal system of L 2 (M) of complex square integrable functions defined 
over M with respect to the measure dpj. 
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It is straightforward to see that the family --™r- n + 2 >— n j s complete in L 2 (M ) with 

respect to the measure d/ij. For any function F € L 2 (M), we define the function / £ L 2 (D) 
induced by the mapping <p in the form / := Foip. The Zernike system is a complete orthonormal 
system in L 2 (D) and then 

oo n 

= E E c ” z ” 

n= 0 m=—n 

a.e., with c™ given in (4). Using that Z™ = K™ o ip we find that 

oo n 

F = E E ■ K « (6) 

n=0 m=—n 

a.e., with 

C = (*C,n^:= [ [ Kn-Fdlij. (7) 

J Jm 


3 Zernike-like systems in polygons 

In this section we particularize the above general formulation to polygonal domains, that is, the 
set M is a polygon of p > 3 sides and radius Rq centered at the origin. We define a := ir/p and 
assume that one of the radius of the polygon is located at an angle a from the positive X axis, 
as it is indicated in right picture of Figure 4 for the particular case p = 8. 



Figure 3: Right picture: octagon of p = 8 sides centered at the origin. One of the sides is 
orthogonal to the X axis and is divided in two equal segments by the positive X axis. Left 
picture: division of the unit disk into p identical sectors. 


The key point here is just to write down a convenient mapping ip from the disk D to the 
polygon M. To this end, we first find the mapping ipo from the sector of the disk S a : = 
{(pcoscj), psin(j))-, —a < (j> < a, 0 < p < 1} C D into the triangular portion of the polygon M 
located over the positive X axis; the triangle of vertexes (0,0) and {Rq cos a, ±i?o sin a), (ipo 
maps the shaded sector of the disk D into the shaded triangle of the octogonon M in Figure 4.) 
This mapping is given by 


(x,y) 


(r cos 8, r sin 9) = ipo(p,cj)) 


Ropcos a 



sin (j> 
cos 4> 


0 < p < 1, —a < 4> < a, 
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where, for convenience, we have used polar coordinates (p, <j)) in the disk. Observe that 6 = <j) 
and r = R^p cos a/ cos 4>. 

Now, in order to find the complete mapping ip : D —>• M, we just need to consider a division 
of the disk D into p sectors of amplitude 2a and repeat the mapping po from everyone of these 
p sectors to everyone of the p triangles that compose the polygon M (see Figure 4). To this end 
we define the function 


U a (<j>) :=4>- 


4> + a 
2 a 


2a, 


0 < 4 > < 2 t r. 


This function is piece-wise linear with slope 1, discontinuous at 8 = {2m + 1 )a, for integer m 
and bounded by ±ct. We also define the variable radius 


Rol{<!>) : = 


i?o cos a 
cos {U a {4 >))' 


We find that i?ocosa < R a ((j)) < Ro■ The function R a (4 >) is an increasing function for cj) £ 
( 2 ma, (2m + l)a) and decreasing for <f> € (( 2 m + 1 )a, 2 (m + l)a), m = 0 , 1 , 2 ,.... 

With these preparations, we can finally define the mapping p that transforms the unit disk 
D into the polygon M: 


(x,y) = p(u,v) = (ui? a (arctan(u/u)), ui? a (arctan(u/u))), 


u 2 +v 2 < 1 . 


Or, in polar coordinates: 

r (p, 4>) = pR a (^), 0(p, 4>) = (t>, 

The jacobian of the inverse transformation 


(u,v)=<p (x,y) = 


x 


0 < p < 1, 0 < 4> < 2ir. 


y 


k ii Q (arctan(y/x)) ’ i? Q (arctan(y/a;) J 
is, in polar coordinates, J a (8) = [R a (9)]~ 2 . From Section 2 we find that the family {”3 r '' +2 ’---> n ) 

x y 


K™(x, y) = Z r n 


R a (axct&n(y/x)) ’ R a (&rct&n(y/x)) 


or, in polar coordinates, 


K™(r,0) = Z ; 


Ra{8) 


,6 1 , 


( 8 ) 


(9) 


is a complete orthonormal system in L 2 (M) with respect to the measure (in cartesian and polar 
coordinates): 

dxdy rdrdQ 

dpj = —fyr, --—7—rwyi d PJ = 


7 ri ? 2 (arctan(y/x)) 


irRKey 


( 10 ) 


More precisely: 


N nR^ (t^) COS (m6), 

-Ny An " 1 sin (m9), m < 0 , 


Kn(r,8) : = 


m > 0 , 


( 11 ) 
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where the polynomials R^ are given in ( 2 ). Also, the family (iZ ™}™ =0 ”3 n+A '" ,n ; with 


m=—n,— n+2,...,n 


*C(*,i/) = 




Z2 a (arctan(y/a:)) n \R a {axctan.{y / x)) ’ Z2 Q (arctan(y/x)) / 


or, in polar coordinates, 


H™(r,6) = 


yrt 

Rate) n \R a (9) 


, 61 , 


( 12 ) 


(13) 


is a complete orthonormal system in L 2 (M) with respect to the measure (in cartesian and polar 
coordinates) 

, dxdy rdrdd 

dp = -, dp = -. 

7 r 7T 

Then, any function F E L 2 (M) can be written, a.e., in the form 


n=0 m=—n 


m 
n 5 


with 




Also, any function F E L 2 (M) can be written, a.e., in the form 


F(x,y) = '£ Y, <? H n, 


n= 0m=—n 


with 

C = - / / Hn( x iy) F ( x >y) dxd y- 
71 J Jm 

Several representative examples, corresponding to various Zernike wavefront aberrations: 
tilt, defocus, astigmatism, coma, trefoil and spherical aberration are represented in Figure 5 for 
the particular case of the octagon, a = 7 r /8 and Rq = 1 (Only positive values of m are shown 
since m < 0 are rotated versions of the same aberration modes). 


4 Zernike-like systems in an ensemble of hexagonal facets 

In this section we consider a set A composed by the union of N equal polygons of p sides. A 
typical example frequently used in the primary phase of SMTs telescopes consist of 18 hexagons 
Mfc, k = 1,2,..., 18 disposed as in Figure 6 . 

Denote by {(a:*;, yk)}k=i, 2 ,...,N the cartesian coordinates of the centers of every one of the N 
polygons Mfc of A. Then, from the previous section, we know that each one of the N families of 
functions 


[iCp)(x, 2 /), k = 1,2,... ,N, n E N U {0}, 


m = —n, —n + 2 ,..., n, 


defined by 


\ K n] (k X x ,y) ■= K n( x ~ x k,y- 2/fc), 










Figure 4: Level sets for the first few polygonal Zernike functions K™(x,y) in an octagon, that 
determine the first classical aberrations. 


with K™(x,y) given in (8) (or (9) in polar coordinates), is an orthonormal system in L 2 (M k ) 
with respect to the measure dyj given in (10) with (x, y) replaced by (x — x k , y — y k )- 

On the other hand, any function F in the set L 2 (A ) of complex square integrable functions 
defined in A may be decomposed in the sum of N functions F k € L 2 (M k ), each one supported 
in each one of the polygons M k : 

N N 

F(x,y) = F(x,y)^2xM k (x,y) = ^F fc (x,y), (14) 

k=1 k=1 

where XM k (x,y) is the characteristic function of the polygon Mf.. From the previous section we 
know that, for every k = 1,2, 3,..., N, 

oo n 

ft'EEBPra 1 * 1 , (is) 

n=0m=—n 


a.e. in Mf., with 




1 

7T 


K™(x — Xk,y — yk)F k (x, y) 


dxdy 


M k 


i? 2 (arctan((y - y k )/(x - x k ))' 


(16) 


Then we find that 

N oo n 

f(i,i) = EE E (17) 

fc=l n,=0 m=—n 

a.e. in A. In other words, any complex square integrable function in A may be written, a.e., 
as a linear combination of elements of the orthonormal system {[LT™]( fc )(x, y)xM k (x, y)} with 
k = 1,N, n = 0,...,oo and m = — n, —n + 2,... , n. 
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Figure 5: Left picture: disposition of the 18 hexagonal facets designed for the primary mirror of 
the James Webb space telescope projected by NASA (picture taken from http://www.quantum- 
rd.com/2010/05/el-james-webb-space-telescope-jwst.html). Right picture: disposition of the 18 
hexagons that conform the set A with N = 18 and p = 6. 


In the particular case considered in Figure 6 of an ensemble of IV = 18 hexagons (p = 6) we 
have that a = 7r/6 and the centers (xk,y k ) of each one of the 18 hexagons M k is located at the 
points 


{x k ,y k ) = V3Ro ^cos ( ^ k ,sin , 


1 < k < 6 , 


(x2k+l,V2k+l) 



3<k<8, 


(x 2 k,V 2 k) = 3Ro 





4<k<9. 


In the appendix we give details about a second example of disposition of hexagonal facets 
for the primary phase of biger SMT telescope considered in the simulations that we show in the 
next section. 


5 Atmospheric turbulence simulation and results 

In this section we test the expansion (17) for a realistic simulation of a wavefront distorted 
by an atmospheric turbulence at the pupil plane of a very large ground based telescope. It 
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is well-known that the random distortions associated to an atmospheric turbulence follow the 
Kolmogoroff statistical model [27]. Roddier [27] proposed an algorithm to generate realizations 
of these distorted wavefronts in terms of Zernike coefficients (or amplitudes). Each realization 
is obtained by generating normally-distributed random numbers with zero mean and a given 
covariance matrix according to the Kolmogoroff statistics. We consider the covariance matrix, 
given in [23]. This matrix, proposed by Noll from energy considerations, is related to (L>/ro) 5 / 3 , 
an overall scaling factor that depends on the telescope’s diameter D and the Fried’s parameter 
ro, which represents the coherence length of the wavefront, or the maximum diameter of an ideal 
diffraction-limited telescope having the same optical resolution (given by the seing disk). 

For this simulation, we consider the pupil of the Gran Telescopio Canarias (GTC), with a 
primary segmented mirror that has a diameter of 10.4 metres. It is located in La Palma (2.326 
metres of altitude) in the Canary Islands of Spain. The primary mirror consists of 36 hexagonal 
facets (disposed as indicated in the appendix). The segmented primary mirror is adaptive, which 
means that both the shape, and relative position of each facet can be changed to compensate 
atmospheric turbulences, and other possible deformations of the telescope structure. For this 
reason, it is crucial to determinate, as much accurately as possible, the wavefront in these 
adaptive optical telescopes. 

The simulated random wavefront follows the Kolmogoroff turbulence model. Using the 
covariance matrix given in [23], with D = 10.4m, and fried’s parameter ro = lm (this is a 
typical value for a reasonably good seing), we have that, in polar coordinates, the fourth order 
approximation to one particular realization of the random atmospheric wavefront at a certain 
instant of time is the function, in phase radians, 

4 n 

/W)T E ( 18 > 

n=l m=—n 

where Z™ = Z™(r,0) are given in (1), the coordinates (r, 0) are the global polar coordinates 
corresponding to the reference system given in Figure 7 (left) and Figure 11: x = rcosO, 
y = rsinO, and the coefficients a™ are computed using the covariance matrix given in [23]. 

We approximate the input wavefront (18) by the expansion (17) up to 4th order: 

36 4 n 

f(r, 6) ~ / approx (r, 0) := EE E (19) 

k= 1 n=0 m=—n 

where the Fourier coefficients [c™]( fc l are given in (16), and (r^, Ok) are the local polar coordinates 
given in the hexagon k: x — Xk = cos Ok , y — Vk = He sin Ok■ 

The function (18) represents an incident wavefront before it crosses the telescope’s pupil. The 
right hand side of (19) is an approximation to the wavefront immediately after passing through 
the telescope pupil. We compute the coefficients [c™]( fc ) in (19) in two different ways. On the 
one hand we obtain [c™y- k > just from their definition: computing the integral (16). On the other 
hand, we use the least square approximation method by fitting the values of the wavefront (18) 
in a discrete data set given by a uniform hexagonal mesh in the mirror containing 25 points in 
each hexagon, as detailed in Figure 7 (left). 

Therefore, depending on which of the two methods we use to compute the coefficients [c™]( fc ), 
we get two different approximations of the function (18) immediately after passing through the 
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Figure 6: Left picture: uniform mesh of the mirror. Right picture: density plot of the wavefront 
(18) on the GTC mirror telescope for the particular realization [27] a™ = {4.25, -9.28, -0.024, 
-0.11, -0.37, 0.28, -0.75, 2.002, 0.49, 0.35, 0.022, -0.009, 0.34, -0.02}. 


telescope pupil: a “least square approximation” and “an integral approximation”. The least 
square approximation allows us to compute approximately the coefficients when we do 

not know the functional expression of the wavefront, but we only know its values at the mesh 
points. The integral approximation gives us the exact value of the coefficients [c™]^; however, 
the calculus of the integral (16) is complicated and requires the knowledge of the functional 
expression of the wavefront. 

The accuracy of both approximations is illustrated in Figures 8 and 9 for the wave front (18) 
with the values of the coefficients a n given in the caption of Figure 7. Figure 8 shows the density 
plots of the exact function (18) and the approximation j a PP rox ( r) given in the right hand side 
of (19) with coefficients computed by both methods. We define the relative error of the 

approximation in the form: 


/ re M) 


/M) - / approx M) 

f(r,e) 


( 20 ) 


where / approx (rj, 0j) represents either, the least square, or the integral approximation. Because 
of the small values of f re (r, 8), in order to visualize graphically the behavior of this function, we 
show in Figure 9, its cubic root. We define the relative root mean squared fitting error in the 
form: 


rrmse approx : = 


M 


1/2 




i=l 


M = total number of points of the mesh. 


( 21 ) 


In Figure 10 we show rrmse approx for some simulations using both methods. 

These results suggest that the accuracy of the estimations is reasonable. The two methods 
used to compute the coefficients the least square and the integral method, seem to be 

equivalent in practice. 
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Figure 7: Left picture: function (18) for the particular coefficients given in the caption of Figure 
7. Middle picture: the least square approximation. Right picture: the integral approximation. 



Figure 8: Cubic root of the relative error function / re (r, #) given in (21), for the least square 
approximation (left), and the integral approximation (right). 



Figure 9: Relative root mean squared fitting errors (21), obtained using the least square and 
integral approximation for 7 different realizations of (18) (7 different sets of coefficients a™). 
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6 Final Remarks 


In this paper we present a rigorous formulation for the general method introduced in [20] to 
obtain complete orthonormal Zernike-like systems in any set M of the plane piece-wise diffeo- 
morphic to the unit disk. It consists of applying a diffeomorphism to the Zernike system on the 
unit disk to transform the circle basis functions into the new system on the desired geometry. 
We then particularize the method to polygons and ensembles of polygonal facets, that are es¬ 
pecially important in the wave theory of aberrations and optical image formation in segmented 
mirror telescopes. In Section 5 we perform a numerical experiment that shows the efficiency of 
the Zernike-like system obtained in Section 4 for the SMTs. 

We believe that the mapping method proposed here, implemented as a change of variables, 
overcomes most of the difficulties and drawbacks of previous methods, and provides a common 
framework, especially well-suited for a unified theory of aberrations. On the other hand, the 
precise physical meaning of the new zernike-like basis K™(x,y) deserves a further investigation. 
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8 Appendix 


A second important example frequently used in the primary phase of SMTs telescopes consist 
of 36 hexagons Mk, k = 1,2,..., 36 disposed as in Figure 11. 

The centers (xk,Vk) of each one of the 36 hexagons Mk is located at the points 


(x k ,yk) = v^3i?o ^cos ’ sin ( 


(k — 1)7T 


, 1 < k < 6, 


y))> 7<fc<12, 


(x 2 k- 7, V2k-7 ) = 2v / 3i? 0 f sin f - + — j , - cos ( - + 

/ / 77 klT \ (77 /c77 \ 

(■£ 2 fc— 18 ; 2 / 2 fc—is) = 3i ? 0 ( cos ( - - — ) , - sin ( - - — j ) , 13 < k < 18, 


( x 3k-38,V3k-3s) ~ 3 V 3 .Ro ^COS 

/— ( ({k ~ 25)7t 

{X 3 k- 55 ,y 3 k- 55 > = V21.Ro I COS I - - - 

( 

/— ( ((k — 31)vr 

\X3k—72 j U3k—72 ) = V21 R 0 I COS I - - - 


6 3 ) 

(k — 19)7t 


, Sill 


(k — 19)7t 


, 19 < k < 24, 


V 3 \\ 

+ arctan I -g- I I , sin 


25 < k < 30, 


+ arctan 
31 < k < 36. 


V3 


, Sill 


(k — 25)ir 


(k — 31)7t 


[V 3 

+ arctan —— 


+ arctan 


V3 N 
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Figure 10: Disposition of the 36 hexagons that conform the set A with N = 36 and p = 6. 
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